Temperature-dependent electronic structure and ferromagnetism in the d = oo 

Hubbard model 
studied by a modified perturbation theory 
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The infinite-dimensional Hubbard model is studied by 
means of a modified perturbation theory. The approach re- 
duces to the iterative perturbation theory for weak coupling. 
It is exact in the atomic limit and correctly reproduces the dis- 
persions and the weights of the Hubbard bands in the strong- 
coupling regime for arbitrary fillings. Results are presented 
for the hyper-cubic and an fcc-type lattice. For the latter 
we find ferromagnetic solutions. The filling-dependent Curie 
temperature is compared with the results of a recent Quan- 
tum Monte Carlo study. 
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Correlations among itinerant electrons are responsible 
for various interesting phenomena like spontaneous mag- 
netic order, the metal-insulator (Mott) transition and 
high-temperature super-conductivity. One of the sim- 
plest but non-trivial models describing -inorrclatcd elec- 
trons on a lattice is the Hubbard modelEl. Studying the 
Hubbard model in the limit of infinite spatial dimen- 
sionsfla d is of special importanae for the construction 
of a dynamical mean-field theoryu. For d = oo the Hub- 
bard model is considerably simplified but nevertheless 
remains non-trivial: It becomes equivalent j-tci an effec- 
tive impurity problem and can be mappecD'a onto the 
single-impurity Anderson model (SIAM), for example. 
For the latter numerically exact solutions can be obtained 
by Quantum Monte Carlo (QMC) calculations. On the 
other hand, it may be helpful to have an analytical (but 
approximative) expression for the self-energy at one's dis- 
posal, which recovers the exact QMC results as reliable 
as possible. This allows a direct calculation of dynamic 
quantities on the real energy axis for finite temperatures 
as well as for T = 0. Furthermore, valuable hints for 
the approximative solution of more complicated lattice 
models not accessible to QMC may be obtained. 

Eog the symmetric case the iterative perturbation the- 
oryQ'B (IPT) is known to give a rather realistic descrip- 
tion of the Mott transition. The IPT employs the self- 
consistent mapping onto the SIAM which is solved by 
means of second-order perturbation theory around the 
Hartree-Fock solution (SOPT-HF). To extend the IPT 
to non-symmetric cases, Kajueter and KotliarQ proposed 



an interpolating expression for the self-energy which for 
arbitrary band-fillings reproduces the atomic as well as 
the weak-coupling limit. 

For the study of ferromagnetism special attention has 
to be paid to the strong-coupling regime. Fortunately, 
the 1/U perturbation theory of Harris and LangelS pro- 
vides rigorous results for U —^ oo: In the first non-trivial 
order beyond the atomic limit, the average dispersions as 
well as the weights-pf the two dominating Hubbard bands 
are known exactlya. One possibility to account for these 
strong-coupling results in the construction of an analyt- 
ical expression for the self-energy, is to ensure that the 
moments 
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of the resulting spectral density A\s^„{E) are correct up 
to n — 3. These moments can be calculated exactly 
by M^a^ = ([/:"ck.,4j+), where C = [-,?{].. At 
each k point of the Brillouin zone the first four moments 
n — 0, . . . , 3 provide four pieces of information that in 
the strong-coupling regime determine the dispersions and 
the weights of the two Hubbard bands. This moment 
approach has been successfully employed beforehand to 
improve upon the Hubbard-I solutiona. ^ 

The approach of Kajueter and KotliarQ reproduces the 
correct moments up to n = 2. The n — 3 moment, how- 
ever, may be of particular importance in the context of 
ferromagnetism since it involves a higher-order correla- 
tion function ( "band shift" ) , 
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the (possible) spin-dependence of which is known to 
favour magnetic order and thereby decisively influences 
the magnetic phase diagramcl. 

Based on these considerations the present authors have 
proposed an iraorovement of the Kajueter-Kotliar ap- 
proach recently^ which is correct up to n = 3 (hereafter 
referred to as "modified perturbation theory", MPT). 
The main purpose of this paper is to study the effects 
of the band shift B„ in the paramagnetic as well as in 
the ferromagnetic phase. We compare our results with 
corresponding findings of previous QMC studies. 

Let us briefly recall the essentials of the theory. Details 
can be found in Ref. QiJ. In the d — oo Hubbard model the 



self-energy is a local quantityaH. Via the mapping onto 
the effective impurity problem, it can be regarded as the 
self-energy of the SIAM provided that the hybridisation 
function of the SIAM, is suitably chosen. According to 
Kajueter and KotliarQ, we consider the following ansatz 
for the self-energy of the SIAM: 
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where E^ {E) is the second-order [U ) contribution 
(SOC) to the self-energy within the SOPT-HF, while a^ 
and b„ are free parameters which will be chosen such that 
the moments (tl) 
an and h„ to bet 



ire correct up to n = 3. This requires 
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ric denotes the occupancy of the impurity level, Ua is 
its Hartree-Fock value, and Jl^ is a fictitious chemical po- 
tential that appears in the definition of the Hartree-Fock 
Green's function. The parameter Jig- is fixed by imposing 
the condition njj — n„ (for further discussions on this 
point see Ref. |o|). Finally, B^ is the Hartree-Fock 
value of the band shift B^^. All expectation values in the 
theory can be expressed in terms of the spectcal density 
and thus can be determined self-consistentlyt2l. 

The modified perturbation theory (MPT) reduces to 
the IPT for small U and is exact in the atomic limit. 
By construction it is fully consistent with the rigorous 
results of Harris and Lange for C/ ^ oo and yields the 
correct moments of the spectral density up to n = S-for 
arbitrary U. The approach of Kajueter and KotliarlZl is 
recovered if B„ and Ba are set to their atomic limit 
values (i. e. B^r — Bo- — Ta = 0). Note, however, that 
in Ref. ^ the fictitious chemical potential Jlu has been 
used to enforce the Luttinger theorem. This implies that 
the theory is intrinsically limited to T = (sec Ref. 10). 

The more unproblematic condition ria = n^ which is 
chosen here to fix Jla- has to be preferred since it allows to 
perform finite-temperature calculations, too. The results 
for r = and near half-filling show that this choice is 
still consistent with the Luttinger theorem (see below). 

Let us first consider the paramagnetic Hubbard model 
on the hyper-cubic lattice away from half-fiUins (n < 1). 
The bloch density of states (BDOS) is given byi poiE) = 
exp(— £'^)/-y/7r. Fig. 1 shows the density of states (DOS) 
for [/ = 4, kgT = 0.138 and different band-fillings ob- 
tained within the MPT. The copparison with the QMC 
results of Jarrell and PruschketJ shows that the MPT 



qualitatively yields the correct results as concerns the 
shifts as well as the changes in height and width of the 
high-energy charge-excitation peaks in the spectrum with 
varying filling. Even for temperatures T ^ the Kondo- 
type resonance shows up at i? « //. This low-energy fea- 
ture is more pronounced in the QMC spectra compared 
with the MPT. We then investigated the effect of taking 
into account the n = 3 moment: For any filling consid- 
ered, it turns out that there are only minor changes in 
the spectra (which would hardly be visible on the scale 
in fig. 1) when setting B^ = B^ . For the paramagnet 
we conclude that the effects introduced by the band shift 
are rather unimportant. 
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FIG. 1. Densities of states for (7 = 4, ksT = 0.138 and 
different fillings on the he lattice. Results of the modified per- 
turbation theory (MPT) in comparison with Quastum Monte 
Carlo (QMC) results from Jarrell and PruschkeEil. 



Fig. 2 shows the MPT-spectra for [/ = 4,n = 0.94 
and different temperatures. We have also calculated 
the uniform static susceptibility in the paramagnetic 
phase by applying an infinitesimally small field: x — 
-§jj(n^ ""■i)|if_+o' ^^ '^^^ ^® s&ev^ in fig. the Kondo 
peak starts growing when the susceptibility differs from 
the Curie-law. The local effective magnetic moment 
(/i^fj (X Tx{T)) is quenched as T ^ 0. This effect 
has been termed "collective single-band Kondo-effect"ll3. 
Furthermore one can see in fig. by comparing the T = 0- 
spectrum with the free ([/ = 0) spectrum that the DOS is 
unrenormalised at the Fermi edge, Paip) ~ Po(mo), which 



is an equivalent foKaaulation of the Luttinger theorem in 
infinite diniensionstS. Let us stress once more that con- 
trary to Ref. |j this has not been enforced via the choice 
for the parameter Jicr- For fiUings well below half- filling 
Pa{l^) is found to be slightly lower than po{iiq). We did 
not find ferromagnetic solutions on the he lattice. The 
susceptibility never diverged for [/ < 6. This is in agree- 
ment with the results in Ref. y_l[ For [/ = oo, however, 
a region of non-vanishing ferromagnetic polarisation is 
obtained within the non-crossing approximationo. 



while at high temperatures we still have Curie- Weiss be- 
haviour. Finally, for n > 0.69, x~^ remains positive for 
all temperatures. This scenario is well known from previ- 
ous studies of the Hubbard model p,conservinglJ as well 
as non-conserving approximationsES. For n > 0.67 the 
true Curie temperature Tc (with Ti < Tc < T2) can be 
found e. g. by means of the Maxwell construction consid- 
ering the H-m isotherms at different temperatures. On 
the other hand, we cannot exclude that the first-order 
transition is an artefact of our approximation. 




FIG. 2. Densities of states iox U = A, n = 0.94 and different 
temperatures. Thin dashed line: [/ = 0. Inset: uniform static 
susceptibility x vs temperature. 



On an fcc-type lattice ferroroagnetism has been found 
recently within QMC by Ulmkeli3. The high weight in the 

BDOsEl po{E) = exp(-(l + y/2E)/2)/ ^ t:{1 + V2E) at 

the lower band edge favours ferromagnetic ordcrta. Sta- 
ble ferromagnetic solutions are likewise predicted by the 
MPT in a wide region of the phase diagram. Fig. 3 shows 
the magnetisation m — n^—rii for [/ = 4 as a function of 
temperature for different band-fillings. For n < 0.67 we 
observe second-order phase transitions. The magnetisa- 
tion curves are continuous and terminate at the respec- 
tive Curie temperature Tc- The inverse static suscep- 
tibility obeys a Curie- Weiss law for high temperatures. 
For lower temperatures with T > Tc a slight curvature 
(9^X~^/9T^ > 0) is observed (not visible in fig. 3). x~^ 
vanishes at T = Tc- The critical exponent of the suscep- 
tibility turns out to be 7 « 1 as expected. We encoun- 
tered numerical difficulties to obtain an accurate value 
for 7. 

For n > 0.67 the phase transitions are of first order. A 
finite magnetisation is still found for temperatures above 
the zero Ti of x~^ (see fig. 3 for n = 0.67). At a temper- 
ature T = T2 > Ti the magnetisation non-continuously 
drops to zero. Between Ti and T2 there is a second fer- 
romagnetic solution with a magnetisation (not shown) 
that vanishes at Ti and that coalesces with the plotted 
magnetisation curve of the first solution at T2. Decreas- 
ing n also results in a stronger curvature of x~^ near Ti 
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FIG. 3. Thick lines: magnetisation for the fcc-lattice, 
17 = 4 and n = 0.2, 0.3, 0.4, 0.5, 0.6, 0.65, 0.67, 0.7, 0.8, 0.85 vs 
temperature. Thin lines: corresponding inverse static suscep- 
tibility x~^ for n = 0.2, . . . , 0.65. Dashed line: n = 0.67. 



The filling-dependence of Tc (for n > 0.67: Ti) can 
be compared with the available numerically exact QMC 
results from Ref. 15 where the Curie temperature has 
been obtained by linear extrapolation of the inverse sus- 
ceptibility. As can be seen in fig. 4, the MPT result 
reasonably agrees with the QMC data. In particular, 
there is remarkable agreement with respect to the maxi- 
mum Tc and the corresponding filling. Contrary to the 
paramagnetic phase, there is a strong effect of the (spin- 
dependent) band shift B^ for the ferromagnet. Both, 
the range of existing ferromagnetic solutions as well as 
the values for Ti and Tc , are strongly affected if we set 
Bcr ~ Ba ■ We find that taking into account correctly 
the n = 3 moment tends to improve the agreement with 
the QMC results (see fig. 4). 

There are no hints for non- Fermi-liquid behaviour in all 
our calculations. The imaginary part of the self-energy 
vanishes quadratically aX E — p in the para- as well as 
in the ferromagnetic solutions for T = 0, even at the 
quantum-critical point fic ~ 0.16 (for U = 4). 

The evolution of the spin-dependent spectral density 
with increasing temperature is shown in fig. 5 for n = 
0.6. Since the self-energy is k-independent, the spectral 
density depends on k via the Bloch dispersion e(k) only. 
For the calculation of the density of states which is also 
shown in fig. 5, the spectral density has to be weighted 
with the BDOS: 



p,{E) - i yde(k) po(e(k))A,(k).(i? - A*) 
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At the temperature ksT — 0.022, where the system is 
nearly fully polarised (m = 0.57), the upper Hubbard 
band in the |-spectrum is missing because there are not 
enough interaction partners. The same spectrum without 
the narrow Kondo-type resonance (at £^ w 0) is expected 
in the non-interacting case. The J,-spectrum shows the 
lower as well as the upper Hubbard band. When increas- 
ing the temperature (second row, ksT = 0.058, m = 
0.32) the upper Hubbard band in the "f-spectrum comes 
into existence. The resonance smears out and starts 
growing in the J,-spectrum. At ksT — ksTc = 0.065 
(third row) where \- and |-spectra coincide, it is still 
visible. Finally, we notice that there is a strong trans- 
fer of spectral weight between the Hubbard bands with 
increasing temperature. In the |-channel the upper one 
gains weight at the cost of the lower, in the |-channel the 
situation is reverse. 
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FIG. 4. Filling dependence of the Curie-temperature Tc 
(T2 for n > 0.67) for [/ = 4 on the fcc-lattice. Solid lines: 
MPT. Dashed lines: MPT results with Ba = bS"^' = 0. 
Thin (solid and dashed) lines: 7i (zeros of X~^)- Ppints with 
error bars: QMC results (zeros of x~^) by UlmketJ. 




FIG. 5. Density plots of the spectral density: energy on 
the y-axis, Bloch dispersion e(k) on the x-axis. f/ = 4, 
n — 0.6 and the fcc-lattice. First row: |-spectral density, 
j,-spectral density and density of states Pcr{E — jj.) (solid 
line I, dashed line: J,) at fcsT = 0.022 (magnetisation 
m = 0.57). Second row: ksT = 0.058, m = 0.32. Third 
row: ksT = ksTc = 0.065. 



In conclusion, ensuring the correctness of the moments 
of the spectral density up to n = 3, is a necessary condi- 
tion to be consistent with the rigorous results of Harris 
and Lange for the strong-coupling regime. In this respect 
the modified perturbation theory not only represents a 
conceptual improvement upon the approach of Kajueter 
and Kotliar, but also yields closer agreement with QMC 
data. The (possible) spin-dependence of the higher-order 
correlation functions {B„) that appear in the n = 3 mo- 
ment is important for ferromagnetic order and has been 
shown to affect critical fillings and temperatures consid- 
erably. On the contrary, there are only minor effects of 
B^ for the paramagnet. Improvement upon the IPT can 
also be expected for the antiferromagnet at half-filling. 
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